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ABSTRACT 

The purpose of this investigation was to determine the feasibility 
of using an analog computer to obtain the time derivative of the tempera- 
ture response of a compact heat exchanger surface subjected to a step 
change in incoming fluid temperature; and to investigate the effect of 
the ratio of flow length to hydraulic diameter (L/D,,) on the heat 
transfer and flow friction characteristics of compact heat exchanger 
surfaces. 

The method of maximum slope developed by Locke and modified by 
Howard to include conduction pa Semeter as used to determine the heat 
transfer information included herein. 

The results show that an analog computer can be a useful tool to 
aid in the collection and reduction of data. Results in the L/D, 


investigation were generally inconclusive and bear further investigation. 
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NOMENCLATURE 


English Letter Symbols 


A Matrix total heat transfer area sq ft 
A. Matrix minimum free flow area sq ft 
A. Matrix total frontal area sq ft 
A Matrix solid cross-sectional area sq ft 
: available for thermal conduction 
a. Fin thickness Ge 
b Flow passage perimeter ( A.) ft 
c Jet contraction-area ratio for dimensionless 
Cc 
circular tube 
Ce Fluid stream thermal capacity Btu/(hr deg F) 
(mc_) 
P 
C. Matrix thermal capacity (Woe) Btu/deg F 
c Fluid specific neat at constant Btu/(lbm deg F) 
P pressure 
C. Matrix material specific heat Btu/(1lbm deg F) 
Dis Flow passage hydraulic diameter ft 
(4r, ) 
E Friction power per unit area hp/sq ft 
c Flow stream mass velocity (m/A_) lbm/(hr sq ft) 
e Proportionality factor in Newton's 32.2(1lbm ft)/(1bf Bec") 
Second Law 
h Surface heat transfer coefficient Btu/(hr sq ft deg F) 
for convection; heat transfer 
power per unit area per degree 
temperature difference 
K. Contraction loss coefficient for dimensionless 
entrance to heat exchanger 
K Momentum velocity-distribution dimensionless 
d coefficient 
K, Expansion loss coefficient for dimensionless 


heat exchanger exit 


k Fluid thermal conductivity Btu/(hr sq ft deg F/ft) 
kK. Matrix thermal conductivity Btu/(hr sq ft deg F/ft) 
L Total matrix flow length ft 
m Mass flow rate lbm/bhr 
P Pressure lbf/sq ft 
P Matrix porosity (A /A,_) dimensionless 
q Heat transfer rate Btu/hr 
R Gas constant (53.35 (ft lbf)/(1bm deg R) 
-air) 
Th, Hydraulic radius (AL/A) ft 
t Temperature deg F 
u Flow velocity ft/sec 
= Matrix volume cu ft 
We Matrix mass lbm 
X Distance along flow passage ft 
from the matrix inlet 
Greek Letter Symbols 
SB Compactness (A/V) sq ft/cu ft 
B Ratio of orifice diameter to dimensionless 
pipe diameter (d /d) 
XQ Difference or Change (time, 
temperature, distance, etc.) 
og Time sec, hr 
i Fluid viscosity lbm/ hr ft 
j~ Density lbm/ cu ft 
Subscripts 
atm Local atmosphere 
ave Average 
£ Fluid (gas, air) 
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Intecél, imbet 

Matrix, mean 

At orifice 

Solid (Matrix material), static 

Standard (temperature and pressure) 

Local conditions 

Inlet conditions (upstream of matrix and heaters) 
Inlet conditions at matrix entrance 


Exit conditions at matrix outlet 


Dimensionless Groupings 


i 


Fanning friction factor; ratio of wall shear stress to fluid 
dynamic head. Plotted as a function of Reynolds No. to geGine 
the surface friction characteristics. 


Colburn j-factor (N fey. This factor plotted vs. Reynolds 
No. defines the surktice heat transfer characteristics. 


Nusselt Number (hD,,/k) ; a heat transfer modulus 


Prandtl Number GALI a fluid properties modulus 
Reynolds Number (4r, Ghae ), a flow modulus 

Stanton Number (n/c), a heat transfer modulus 
Number of heat transfer units (hA/me |) 


Longitudinal iieat conduction parameter for solid 
material Soon 5 oy 


Time parameter (hAG/W oc.) 
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t Introduction. 

One of the methods used to determine the heat transfer characteris- 
tics of compact heat exchangers, that is, those heat exchangers with a 
heat transfer surface area to volume ratio greater than 200, is the single- 
blow or transient test technique developed by Locke (12) and further modi- 
fied to include the effects of conduction by Howard (6). This method has 
been used by several previous investigators at the Naval Postgraduate 
School to evaluate the heat transfer and flow friction properties of 
several different materials and geometries. 

In the single-blow transient testing technique, a heat exchanger 
matrix is subjected to a step change in fluid temperature. The response 
of the exit fluid temperature is monitored following this step change. 

The maximum rate of change of the exit fluid temperature is uniquely 


related to the dimensionless number of heat transfer units, Ne 


<= hA/me 
Experimentally, this response is recorded as a temperature-time curve and 
the maximum rate of change determined from the maximum slope of this curve. 

Location of the maximum slope of this curve has been determined by 
a "cut and try" approach. This method leaves some uncertainty as to the 
exact location of the maximum slope and it was felt that an analog com- 
puter with a differentiating circuit could be used to either determine 
the value of maximum slope directly or indirectly by locating the inflec- 
tion point on the response curve which would be plotted simultaneously 
with the computer output. Therefore, one of the objectives of this in- 
vestigation was to evaluate the use of the analog computer as a differ- 
entiator operating directly on the output of a thermocouple. 

A second objective was to investigate the effect of the flow length 


to dydraulic diameter ratio, L/D» on the flow friction and heat transfer 


characteristics of a compact heat exchanger surface. 


2. Summary of Theory. 

AS Background 

An analytical solution to the "single blow problem" was first 
presented by Anzelius in 1926. Subsequent work was then done by Nusselt 
in 1927, Hausen in 1927 and 1929 and by Schumann in 1929. (14). Schumann's 
solution to the problem involved an explicit solution of the differential 
equations in the form of two infinite series using Bessel functions. 
These solutions generated a family of theoretical curves which were used 
as a standard. A plot of experimental data was then compared with the 
theoretical solutions, and the theoretical curve best matching the ex- 
perimental data was used to determine the desired information. 

Locke (12), differentiated Schumann's solutions and observed that 
the maximum slope of the exit fluid temperature response curve was a 
unique function of Neu! 

Howard (6), using a digital computer and finite difference techni- 
ques determined a series of solutions for Neu vs. maximum slope that in- 
cluded the effects of longitudinal conduction, which had been assumed 
equal to zero in Schumann's solutions. 

B. Theory. 

The following theory treats the problem of determining, for a 
gas flowing through a porous solid, the temperature of the fluid and the 
solid as a function of position and time after the incoming fluid has 
been subjected to a step change in temperature. 

The following assumptions are made: 

1. Fluid properties are independent of temperature. 


2. Fluid flow is steady. 
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3. Homogeneous porous solid. 

4. Thermal conductivities of the solid and the fluid are in- 
finite in the direction normal to flow. 

5. The thermal conductivity of the fluid is zero in the 
direction of flow. 


The following analysis is based on the energy balance shown below: 


m™epte (> hb -t,)dxC> ep (t, + 2% x ) 


-%A, 2 ——> > ~ he As (Bes Zh) 
bate 


The initial conditions and the boundary conditions for the analysis 
derived from an energy balance on an element of the porous solid are: 
1. The matrix is initially at a uniform temperature. 
2. There is a step change in the incoming fluid temperature 
at time equal to zero. 
3. The matrix boundaries are abiabatic. 
Let: t. = temperature of the solid. 
te = temperature of the fluid. 
The heat rates from the energy balance are: 
1. Energy absorbed by the solid = ~, Asc, (as Si) Clone 
2. Heat transferred to the solid by convection = 454 (t,-,) dx 


3. Heat transferred from the fluid by convection = 


of, 


4. Heat transferred in the solid by conduction = 


-k. As os x 
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The energy balances are: 
Leia = (2-1) 
m cp Sr dy the bya tee ities Nucla Ves 
2. Solid: (2-2) 
pAscs 28s dus hbltz-te)olerks As XE ay 


Let 7 = the generalized time variable (dimensionless) 





T = (hA/Woc.)(@ - Wx/mL) (223) 
where: 
h = unit conductance for convection heat transfer (Btu/hr sq ft deg F) 
A = matrix total heat transfer area, (sq ft) 


Wece = matrix thermal capacity, (Btu/deg F) 


6 = time, (hrs) 


We = mass of fluid in matrix, (1bm) 


mass flow rate of fluid (lbm/hr) 


e 
m 


x = distance along flow passage in direction of flow measured from 
the inlet, (ft) 
L = total matrix flow length, (ft) 


Multiplying the second term in (2-3) by eile, and rearranging leads to 


the equation 


T= hA o aad hAx ( we Cp 
Ws cs mLep Wse;, 
where the second term may be ignored because nee is much less than Wc 


when the working fluid is a gas. Thus the time parameter reduces to 





7 ex fA 4 (2-4) 
We S 


Let z = the generalized position variable (dimensionless) 


z= (hA/mc )(x/L) = N- x/L (2-5) 
p tu 
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and z ® ne at the matrix exit where x = L. 
Let A= the longitudinal conduction parameter (dimensionless) 


AS kA, / ine L) (2-6) 


where; kK, = thermal conductivity of the matrix (Btu/ hr ft deg F) 


Solid matrix cross-sectional area available for thermal 


A 
S 
conduction (sq ft) 


Substituting the above dimensionless groups into the energy balance 


equations (2-1,2) and rearranging, the equations become: 


Fluid: Qtr e Crary (2-7) 
at re 
Solid: ae \ Neu ae + (ee ta) (2-8) 


If thermal conduction is assumed zero in the direction parallel to 


flow, equations (2-7,8) simplify to: 


ot 

Fluid: <tr = f.- -t, (2-9) 
ot 

Solid: Ss =f,-f (2-10) 
ee sls 


Schumann's solutions to these equations are: 


oo 
- (2+7) 
t,.-T; YT ae 
fast) aye) 2 (2) aa 


n=/ 


(2-12) 


~\ 00 
(ts -t;) ie ates mtg n dd” ; 
—) /-=e but Tew” ( L(2ivzr)) 


N=oO 


Locke (12) then differentiated these solutions for constant Bea to 


arrive at the following equation for the slope of the heating curve: 
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d [Sa 2 | 
Gp oie Ney — gale Calne e 


d ay, V Neu 7 


-(N,,+ 7) _— 


The response of the downstream fluid temperature is used at x = L 
where z then equals Neu and ty equals tro: 

At low flow rates the effects of longitudinal conduction cannot be 
neglected and these effects were taken into account by Howard (6) who 
applied a finite difference method to equations (2-7) and (2-8) to obtain 
solutions, with the aid of a digital computer, for Ty es as a function of 
both maximum slope and conduction parameter r . Howard's results were 
plotted and tabulated so that by calculating Xr and measuring the maxi- 
mum slope of the heating curve, the graph or the tables could be entered 
to determine the corresponding values of Neu’ See Figure 1 and Table l. 

C. Differentiator. 

Determination of the maximum slope of the exit temperature has 
been done by a visual procedure by previous investigators of this problem. 
Slope was found by sliding a straight edge along the trace of the tempera- 
ture curve until the point where slope appeared to be a maximum was found. 
The tangent line to the curve was drawn at this point and the value of 
slope calculated. It was felt that this method left some doubt as to 
the accuracy of the value of maximum slope and it was decided to use an 
analog computer with a differentiating circuit to determine maximum slope. 
Ideally, the computer output should give the value of slope desired 
directly, but if this method should turm out to be inaccurate, the peak 
of the derivative curve should allow the point of inflection on the 
temperature response curve to be accurately determined, assuming no tim- 


ing errors exist. 
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When using the analog computer to solve differential equations, the 
computer is normally programmed to solve equations by integration, avoid- 
ing differentiation whenever possible because of the inherent noise pro- 
blem. If noise is represented by a sinusoidal such as Asin (G)t), 
where G) equals the angular velocity, then its derivative becomes Aa 
cos (GJt) and the magnitude of the differentiated noise signal is direct- 
ly proportional to frequency.With the output of the thermocouples used to 
monitor temperature response in the vicinity of five millivolts and less, 
high frequency noise signals could easily obliterate the data signal 
after differentiation. 

Since the signal in the circuit was D.C., and no alternating cur- 
rent was present, it was felt that the noise problem was minimal and an 
analog differentiator could be used. 

Several circuits are available for analog differentiation (5), the 


simplest being the inverse of the integrating circuit: 





¢ 
R 
e 
e. ~ 
OPERATIONAL AMPLIFIER 
Integrating Circuit 
R 
Cc 
coal a P e. 
— go 
ce 


OPERATIONAL AMPLIFIER 
Differentiating Circuit 


By summing currents into the amplifier (current through the amplifier 


assumed zero) the following differential equation describes the circuit: 
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ae, . 
Cc = -+ Se =O or ere a aa 
AT R t Re ~° 


where if we let e. be the output of the downstream thermocouples, then 
the output of the circuit (e, above) is equal to the derivative, modified 
by the coefficient (-RC). Taking the Laplace transform and putting the 
equation into the s domain, where 
s = a complex variable 
E (8) = - RC s E, (s). 
Letting the product of RC equal to one reduces the equation to: 
E fs) = - § E,(s). 
Recalling that the transfer function E (8) /E, (s) = gs is the deriva- 


tive, the equation for the circuit can be rearranged to give: 


ant —RCs oY c= _ 
Ej (s) Ei (s) 


when the RC product equals unity. This circuit is the true derivative 
and is not satisfactory for use if any noise is present, due to noise 
amplification. 

Additional circuits that attempt to overcome the noise problem are 
combinations of differentiators and first order, low pass filters. An 
example of such a circuit is the basic differentiator with additional 


resistors or capacitors, shown below with their corresponding transfer 





function: 
R 
R2 
cee C fa 
eo-Ww— e, e, o—_—__| ©. 
E.G) Rac E Cs) 
—— | . S245. © ara RCs 
E; Cs) RiCsel Eid” RCas-+] 
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A drawback of the one amplifier circuits shown above is the necessity 
of knowing the frequency of the noise and precalculating the required 
resistor and capacitor values based on eliminating this noise. 

The modifications of the basic differentiator are referred to as 
approximate differentiators. Several of the circuits were tried in the 
course of this investigation and were found unsuitable. One circuit that 
was satisfactory and is recommended for general use (5), can be found in 
Jackson (7). This circuit combines an adjustable low pass filter with 


the differentiator and is as follows: 


y=-f/zdt <3 


INTEGRATOR 
—(X#AZ+y) +(xtaa+y) az mw GK 
x i at 
U A 
SUMMER SIGN CHANGER 
az 
Potentiometer 


The differential equation that follows describes the problem to be 
solved, that is, an approximation to the derivative. It is: 


(/- a) aE +2 = aH (2-14) 


If (a) is allowed to approach one, the term containing dz/dt diminishes 
and the equation becomes z = dx/dt when (a) is equal to one. Integrat- 
ing equation (2-14) results in an equation that describes the analog cir- 


cuit shown above: 


Za) 


cE 
(i~-A)z+| ZAT = X (2215) 


Or, rearranging 


t 
Z= X*QA2Z-|/zdt 


oa 


(2-15a) 


Taking the Laplace transformation and putting (2-15a) into the s 


domain results in: 


aS) ie X (Sd) + a 2 (Ss) ~+ Z CY rok; (2-16) 


Zi9(s-a+t) = X(s) (2-16a) 
The transfer function for this circuit is: 


ZiGs) gla OF) A Sh (2-17) 
X (5) Fe (S) (1-Q@)se/ 


which reduces to: 


E. (5) in 
EG (s) 


when (a) is equal to one. The features of this circuit which add to its 


(2-18) 


convenience are the potentiometer, which allows adjustment of the low 
pass filter to filter out the desired frequencies without precalculation 


and the sign changer, which provides the derivative with the proper sign. 


D. Effect of L/D,,- 
There is little information published on the effect of L/D, - 
Analytical solutions exist for simple geometries such as circular tubes, 
annuli and parallel planes for the laminar flow fully developed tempera- 


ture and velocity profiles. For complex geometries however, the 


oe 


differential equations become too difficult to solve. From circular 
tube theory it is known that the hydrodynamic and thermal entry lengths 
are a function of Reynolds No. and x/D, where x is the dimension in the 
direction of flow measured from the tube entrance. Kays and London (9) 
list solutions for fully developed velocity and temperature profiles for 
simple geometries that hold for L/D,, greater than 100. For fully devel- 
oped flow, the ratio of the Stanton No., Prandtl No. product to the fric- 
tion factor is a constant. This means that on a plot of f and j versus 
Ne the two factors should plot as a pair of parallel lines. Since fric- 
tion factor is inversely proportional to Reynolds No. in the laminar 
region, for an arbitrary long tube, friction factor approaches a slope 
of -1.0., and the corresponding j - factor does also. 

Kays (9) has a plot for triangular flow passages where he shows an 
analytical solution for an isosceles triangular shaped passage with L/D, 
at infinity. Here the slopes are equal to minus one. 

By using several cores of the same material and flow passage dimen- 
sions and geometries, with only the length and subsequently L/D,, vary- 
ing, one should be able to experimentally determine the ratio of L/D, 
where slopes of the f and j curves are equal and where the flow becomes 
developed to the point that the exit and entrance effects are negligible 


in comparison with the effects of the internal passage. 
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i Experimental Techniques. 

The experimental apparatus was designed to conform to the idealiza- 
tions stated by Howard (6) so that Howard's conduction parameter curves 
could be used. Briefly, these idealizations are: 

1. Fluid flow in the matrix is steady with a uniform velocity 
and temperature profile at any cross-section. Matrix thermal conductivity 
is infinite normal to the flow direction and finite in the direction of 
flow, making the problem one dimensional in space. These idealizations 
have been met through a specially designed entrance nozzle, flow straight- 
ening screens and the heater wires distributed evenly across the flow 
channel. Conductivity is accounted for in r , the conduction para- 
meter. 

2. Large matrix thermal capacity in comparison with that of the 
contained fluid. Use of a gas as the working fluid satisfies this ideal- 
ization and eliminates time dependent terms from the equation for the 
Fliniade, 

3. Constant and uniform thermal properties of the fluid and 
the matrix. 

4. The convective heat transfer coefficient is some suitable 
average and remains uniform and constant. The idealizations of steps 
three and four are met by restricting the temperature changes to small 
values (about 20 deg F above ambient) where the variation in properties 
can be neglected. 

5. A step change in fluid temperature is imposed at time 
equal to zero, after the matrix and entrained fluid have reached a 
steady state temperature. The heaters are of 0.0031 inch diameter nich- 
rome wire with a response time of less than ten percent of the matrix 


response, giving essentially a step change. Heater time constant for 
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a flow rate of 1,000 lbm/hr is 0.0425 seconds compared to a 1.4 second 
time constant for the matrix response. At low flow rates, about 25 lbm/ 
hr, the heater time constant is about 0.3 seconds which compares to a 
time constant of 57 seconds for the system response. Therefore, the 
heaters provide a very close approximation to a step function. 

A schematic diagram of the experimental system is shown in Figure 3, 
which indicates the location of instrumentation. Figures 4 and 5 are 
photographs of the system. Flow is induced into the system through an 
inlet bell, followed by a flow straightening section, the heaters, the 
matrix test section, and a flow measurement device. Pressure taps mea- 
sure the static pressure at the inlet to the test section, pressure dif- 
ferential across the test section, static pressure in the pipe before the 
flow metering orifice and pressure differential across the orifice. 
Thermocouples are located in the inlet to the apparatus, between the 
heaters and the matrix, immediately downstream of the matrix, and in 
the pipe preceding the orifice. 

Pressure is measured by draft gage for low flow rates and water mano- 
meters for the flow rates over about 350 lbm/hr. Atmospheric pressure is 
measured on a standard mercury barometer. 

The temperature response recorded is the difference in temperature 
between the inlet temperature, th» and the matrix exit temperature, t,: 
A more complete description of equipment is found in Appendix A. 

A test run is accomplished by predetermining desired flow rate, cor- 
responding approximate pressure drops and required heaters. Air is drawn 
through the apparatus, flow is adjusted until the desired pressure drop 
across the orifice is attained and the heaters are energized. The matrix 
and the heated air are allowed to come to steady state temperature, then 


power to the heaters is secured and a recording of (t,-t,) as a function 
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of time is made on one channel of a dual channel strip chart recorder. 
The thermocouple output is fed to a differential amplifier, amplified by 
1000:1, then fed to the differentiator on the analog computer. The dif- 
ferentiated signal is then led to the second channel on the recorder where 
it is simultaneously plotted with the time temperature response curve. 
The following information is recorded for each run: inlet static 
pressure, matrix pressure drop, orifice static pressure, orifice pres- 
sure drop, atmospheric pressure, air temperature at the orifice, diameter 
of the orifice plate, ratio of orifice diameter to pipe diameter, chart 
speed of the recorder and corresponding scale factors for the recorder. 
After completion of all desired runs, the values to compute slope are 
taken from the recording traces and included on the data sheet. The 
data sheet layout conforms to the data input section of a digital comput- 
er program that reduces all data and calculates the desired information. 
The derivative was calibrated in the following manner: First, a 
known potential from a potentiometer was used to calibrate the channel of 
the strip chart recorder used to record the derivative. This calibra- 
tion was made in the range of 100 millivolts to 5 volts; which was the 
range of the output of the differentiator on the analog. Again, using 
a known potential, the amplifier was adjusted to give exactly 1000:1 
amplification. A ramp function was set up on the analog, and the output 
of the ramp recorded by the calibrated channel on the recorder. An ac- 
curate measurement was made of the slope of the ramp function from the 
chart trace. The ramp function output was then fed to the differentiat- 
ing circuit and the result recorded on the recorder. From equation (2-14) 
or (2-17) it can be seen that if the potentiometer in the differentiat- 


ing circuit is set equal to one, the circuit reduces exactly to the 
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derivative. However, this was not the case in actuality. 

It was found that as (a) approached one the circuit became more 
unstable, showing up as large oscillations and high overshoot about the 
value of the derivative, and when equal to one, all stability was lost 
and no value for the derivative could be determined. By using the pre- 
viously measured value of slope for a standard, the potentiometer was 
adjusted as close to one as possible, about 0.95, while holding oscil- 
lation and overshoot to a reasonable amount. This resulted in a final 
value of the derivative equal to 99.2% of the value used as a standard. 

Noise is a problem in any derivative circuit and the circuitry used 
for this investigation was no exception. The heating system used in the 
past provided current to the heaters from an A.C. source with voltage control- 
led by a Variac autotransformer. The thermocouples used in the test ap- 
paratus are non-shielded and no shielding was provided on the leads from 
a common terminal strip to the outlets located on the base of the equip- 
ment. In addition, no system of grounding was used in the thermocouple 
circuit. The Moseley dual channel, strip chart recorder used has the 
circuitry to filter out spurious noise signals when this noise was imposed 
on the output thermocouple signal, but when this signal was differentiat- 
ed, the noise level was great enough to completely mask the derivative. 

It was then discovered that much of the interference stemmed from in- 
duced 60 cycle noise from the heater circuitry which is in close proxi- 
mity to the thermocouple circuitry. This was eliminated by using D.C. 

for the source of power to the heaters. There is still considerable noise 
in the system, but the filter on the Moseley recorder is satisfactory to 
handle most of it. For example, on a run where the derivative measured 
178 millivolts/second, the measured noise was 3 millivolts/second peak 


to peak. 
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Amplification of the thermocouple response was necessary so that 
the derivative trace could be plotted, especially at low flow rates. 

At low flow rates the derivative is often less than 0.1 mv/sec which is 
only 10 percent of the lowest scale on the recorder. In order to boost 
the analog output to the more usable ranges of the recorder, an amplifica- 
tion of 1000:1 was used. 

With the amplifier and the analog computer circuits paralleling the 
direct input to the recorder, and because both circuits used the same 
source voltage from the thermocouples, which was only three millivolts 
it was felt that possibly there may have been some interaction between 
systems. Several comparison runs were made to see if this situation did, 
in fact exist. Data was collected for the same flow rate for the follow- 
ing conditions: 

1. Amplifier and analog computer in the circuit. 

2. All components separated from the circuit and direct 
current used for heating alone. 

3. All components separated from the circuit with alternating 
current used for power to the heaters, which duplicates 
the method of previous experimenters. 

All methods of determining maximum slope yielded values of ae 


within the range of experimental accuracy. 
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4, Description of Test Matrices. 

The heat transfer surfaces used for this project were all of a 
similar geometry and construction, consisting of triangular or corrugat- 
ed fins made from solid sheet stock, and splitter plates made from the 
same material. 

The matrices used during the evaluation of the analog technique had 
been evaluated previously by Ball (2) and Bannon (3). These were used to 
provide a correlation on results and to insure proper experimental techni- 
que was used, that is, that reproducible data was obtained. 

For the investigation of L/D,, effects, a matrix with a small hydrau- 
lic diameter was required so that a large value of L/D, could be obtain- 
ed within the limiting dimensions of the test apparatus. Fortunately 
a crimping roller that produced a forty fin per inch by eighteen mil 
high triangular fin was available. One mil brass shim stock was used to 
construct the fins and splitter plates. Test cores were made to the 
following lengths: 1/2 in, 3/4 in, 1.0 am., 1% ins, 2% im. Jatt 340 im. 
A stainless steel reference matrix which had been constructed with the 
same equipment and had previously been tested was available for use as 
a comparison standard for flow friction results. 

Further information on core geometries and properties is shown in 


Figures 6 through LO. 
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Die Presentation of Results. 

For each matrix tested, the heat transfer and flow friction char- 
acteristics have been computed, plotted and tabulated. Computed results 
are shown tabulated in Tables II through XI, and are plotted as Colburn- 
j vs Reynolds Number, and Fanning friction factor vs. Reynolds Number. 
Reynolds Number was calculated on the basis of hydraulic diameter for 
each matrix. 

For the investigation of L/D, effect, six cores were tested. Fig- 
ures 15 through 20 show individual data from each core and Figures 21 and 
22 are compilations of f and j data respectively from all the cores. 

Three cores tested by previous investigators in this project, Solar 
No. 1,Solar No. 4 and the stainless steel plate-fin reference were used 
for a basis of comparison of the analog technique. A plot of "f" and 


"i" vs. N. for these cores is shown in Figures 11, 12 and 14 where the 


R 
results obtained by the use of the analog are compared with previously 
established information. 

During the investigation it was felt that the response of the 
Moseley recorder was not sufficiently rapid at higher flow rates to pro- 
duce accurate results. A comparison was made between the Moseley record- 
er and a Brush recorder using a simultaneous recording technique for the 
testing of one core sample. The results of this comparison are shown in 
Figure 24. 

One new core was evaluated during the course of this investigation. 
Designated Solar No. 6, this core was a stainless steel, triangular fin, 


plate-fin construction with a small hydraulic diameter (0.00126 ft) and 


a fin height of 0.022 inches. Data for this core is shown in Figure 13. 
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Gr, Discussion of Results. 

For the determination of the practicability of using the analog 
computer, several cores previously investigated were chosen to act as 
standards of comparison. These were Solar No. 1, No. 4 and a stainless 
steel reference matrix. Figures 11, 12 and 14 show f£ and j data versus 
Reynolds Number for both methods of determining maximum slope: direct 
reading of analog computer output and the plotting of the tangent line 
at the point of inflection. Plotted on these figures are the data of 
previous investigators at this facility. It was hoped that this data 
could be duplicated so that an accurate comparison could be made. In 
addition, in the L/D, investigation, the data presented shows the results 
of both methods used. 

Examination of Figures 11, 12 and 14 show that flow friction results 
matched exactly and that the heat transfer data match in the region of 


higher Np» where N. is generally above 75. There is a disagreement in 


R 
the lower flow rate region where the data presented ceases to follow the 
expected straight line behavior, but reaches a maximum and then decreases. 
However, closer examination shows that the information presented by both 
methods of determining maximum slope do agree and in fact, several points 
actually coincide. Inspection of the data traces while reducing data 
showed that the cores with a larger hydraulic Pata than that used in 

the cores for the brass L/D,, investigation, namely Solar No. 1 and No. 4 
indicate two possible maximum slope values when the flow rates are over 
about 500 lbm/hr. The trace of the derivative of the response curve shows 
an initial peaking at time equal to 0+. This curve then goes through a 


minimum point and goes through another peak before it decays to zero. 


If slope is computed at both of these peaks and the corresponding 
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j-factor plotted, it is seen that the values at the later peak follow 
the trend of the earlier data, where no ambiguity exists, while the 
initial peak gives erroneous results. Kohlmayr (11) has extended the 
maximum slope method to include the effects of deviation from the ideal- 
ization of a step input. In his paper he describes the double peaking 
of the derivative of the generalized heating curve, and contends that 
the later peak is the correct value of maximum slope to use. The simul- 
taneous plotting of the derivative gives an immediate indication of the 
location of the exact maximum slope and in this case is an aid in deter- 
mining the proper value of the heat transfer coefficient. 

When flow rates increase, the two peaks of the derivative curve move 
closer together until they are indistinguishable. It was felt that the 
recorder in use had a response time too slow to distinguish between these 
values; in order to check out this contention, a simultaneous recording 
was made on the Brush recorder, which has a more rapid response time 
than the Moseley recorder, 0.0035 seconds to 0.23 seconds respectively 
for rise time. No discernable difference between traces could be noted. 
Actually, at the start of the investigation it was thought that the 
initial peaking was due to a transient response that exceeded the capa- 
city of the recorder to handle, it was felt that the use of the Brush 
recorder would differentiate between signals. It was later in the course 
of this investigation that Kohlmayr's results became known. It is in- 
teresting to note that at the upper flow rates, about 1000 l1bm per hr., 
the thermocouple time constant is 0.43 seconds compared with a time 
constant for the downstream -Bamoiee of 1.4 seconds. This could possi- 
bly be a source of error at high flow rates. 

The investigation of the effect of L/D,, on flow friction and heat 


transfer did not produce as good a result as hoped for. Above a value 


a2 


of L/D,» hopefully to be determined, fully developed flow exiats and 

the plot of j and f data should be two parallel lines with the same slopes. 
In other words, j/f should be a constant. For friction, since f is in- 
versely proportional to Np» the slope of the f curve should approach minus 
one as the limiting value as L/D, becomes arbitrarily large. The slopes 
of the £ and j curves were calculated and compared for each value of L/D,,- 
It can be seen that instead of approaching a limiting value the slopes 
went through a maximum with the 1.5 inch core, then decreased somewhat. 
All cores were made by this experimenter except the 3.0 inch and the 1.5 
inch cores. Each core was constructed of the same number of fin and 
splitter plates, and the dimensions of the individual pieces were con- 
stant except for flow length, which varied according to the size core 
under construction. In spite of these precautions, the overall stacked 
height of each core varied sufficiently to affect the geometric constants. 
With only length changing and all other factors identical, porosity, com- 
pactness, frontal area, free flow area and conduction area should have 
been identical for all cores, but all varied. 

The core material, brass shim stock, was prepared by slicing on a 
paper cutter. This produced burring of the cut edges, which in turn af- 
fected the flow characteristics. The previously constructed cores did 
not have burred edges. 

Brass was used primarily because a quantity was on hand and the addi- 
tional amount needed was easily obtainable from local suppliers. It was 
possible to construct and test cores from material on hand while waiting 
for additional materials to arrive. Unfortunately, brass has a high 
thermal conductivity which gives a’ hich value to conduction parameter at 
low flow rates. In addition, conduction parameter is inversely proportion- 


al to flow length and flow rate. These three constituents combine to give 
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high values of conduction parameter at low flow rates and with low flow 
length cores. An inspection of Howard's solutions (Figure 1, Table I) 
show few solutions in the high » region, that is, A > 0.1. 
With high flow rates, maximum slope decreased and the solutions gave 
values of Nig in the range of 1.0 to 3.0, where the possible error is 
greatest (Figure 2). 

While these problems do not account for the "hump" noticed in the 
j-curve, it does give possible reasons for the spread of data, especial- 
ly for the shorter length cores. In retrospect, it appears that a study 
of L/D, effects could better be done by using a low r, core to obtain 


data for larger L/D,, ratios and a dimensionally similar but with a 


larger r, core to obtain data for small ratios of L/D,,5 preferably 


h 


constructed of a low conductivity material. 
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? Experimental Uncertainties. 

There are several sources of possible error which arise from not 
exactly meeting the idealizations and boundary conditions set forth 
previously. These possible error sources are difficult to accurately 
assess a numerical value to, and will not be discussed further in this 
section. 

Error sources which can be evaluated include uncertainties in 
physical constants, inaccuracies in the determination of geometrical 
constants and inaccuracies in instrumentation. Single runs were usually 
made for each data point except where the validity of a run was in ques- 
tion and a check run was made. This system of taking data does not allow 
a large store of statistical data to be collected where inference to error 
may be made from statistical techniques. In view of these limitations, 
the technique for determination of uncertainties developed by Kline and 
McClintock (10) was used. 

The physical constants for the core materials tested were previously 
qeremnined by Ball (2) and Bannon (3) and the uncertainties listed in 
those references are: 


Physical constants 


k + .5% 
Ss. 

Cua oe 
ee 

c + .5% 
p= 

Np, + 2.0% 


fot 1.0% 


Similarly, uncertainties in dimensions have been determined to be: 
+ 1.0% 
a, Ar. Aa» A. + 1.0% 
L + 0.5% 


W << 0.1% 
s 
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The weight of the matrix was measured on an analytical balance and any 
uncertainty can be considered negligible compared to other physical 
dimensions and constants. 

Inaccuracies in instrumentation were generally established to be of 
the order of one half of the smallest division of the scale on the in- 
struemnt in use. The largest source of inaccuracies were in the mano- 
meters used for pressure measurement. All temperature recordings for 
conditions at the matrix were differences in temperature and were record- 
ed in inches on the strip chart recorder. The measurement of the tempera- 
ture at the orifice was obtained from the reading of the thermocouple 
output on a Rubicon potentiometer. The estimated error of the potentio- 
meter was estimated as approximately 0.1 deg F which reduces to about 
0.2% uncertainty. 

The range of pressure measurements was such that several different 
measuring instruments were required, from draft gages at low flow rates 
to a 120 inch water manometer at high flow rates. Because of this range, 
the maximum uncertainty among all readings was used for the analysis. 
These uncertainties are listed below: 
ea 12257 
A P, + 1.25% 

Ap. + 1.70% 
a 
? et 0.005"Hg (negligible) 

Uncertainty in measuring maximum slope is estimated to be 2.0% (14). 
This uncertainty in maximum slope is then considered with Neu and Xr 
to determine the error in Neu from Figure 2. Close examination of Fig- 


ure 2 shows a large uncertainty in Neu for large values of conduction 


parameter at high Neu and for all values of conduction parameter when 


36 


N is near 2.0. 
tu 


Two examples of uncertainty in Colburn j-factor have been cal- 


culated, one for high N. and one for low Ne (14). Using Figure 2 and 


R 


the methods of reference (10), the following uncertainties in Colburn j 


have been determined: 


me SiO Xr = 0 _ error = + 7.2% j error 7.5% 


i 
I+ 


It 
4. 


Ness 25.05 A 008 Ney error = + 10.0% j error 10.2% 


By similar analysis, uncertainties for nm, Ne and £ were determined 


to be 1.0%, 2.3% and 4.3% respectively. 
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On Conclusions. 

The use of a differentiating circuit on an analog computer can be 
used to obtain the maximum slope of the generalized heating curve of the 
transient testing technique for evaluating heat transfer properties of 
compact heat exchanger surfaces. The computer output can be used direct- 
ly to determine maximum slope and indirectly to locate the correct point 
of inflection when multiple maximum slopes occur. 

The effects of L/D, bear further investigation. The results obtain- 
ed in this investigation were inconclusive and at best represent heat 
transfer and friction characteristics of the individual cores tested. 
Better control of geometry and construction is needed to eliminate error 


caused by geometrical inconsistencies. 
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S$. Recommendations for Further Study. 


An investigation should be made in the low N. region to determine 


R 
whether the maximum in the j-data is actually there or if an error in 
technique exists. 

Further investigation of L/D, effects can be made. It would be 
recommended to obtain suitable core samples from an established manu- 
facturer to insure conformity between samples. For investigations in 
the low L/D, range, a core with a similar geometry but larger rh should 


be obtained as the small Th brass cores used were too easily affected 
by high flow rates. 

The cyclic test technique developed by Traister (17), should be 
used for low L/D, cores where Ne is in the maximum error range. 

Further extension of Howard's conduction parameter curves in the 
high New? high lambda region is needed. Insufficient data exists in the 
region where lambda is between 0.1 and 10.0 and is at high Neu’ The 
subroutine which interpolates in this region in the data reduction com- 
puter program often cannot handle these values and gives incorrect re- 
sults. 

A system of shielding and grounding should be accomplished along 
with installation of thermocouples which have a response time less than 
that of the Nichrome heaters. This would make the heaters the control- 
ling factor in the response to a step function rather than the thermo- 
couples which at present cannot handle higher flow rates accurately. 

It is suspected that eHE attenuation noted by Traister in the cyclic 
technique at high flow rates and increased frequency may be due in part 
to the slow response of the thermocouples. Shielding the thermocouples 


and installing a ground system would help alleviate the electronic noise, 
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especially that produced by the linear accelerator located directly 


underneath the test apparatus. 
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eight, W (1bm) 0.4956 


Hicure 10. Physical and Geometric Froperties of Brass Core L=1.5" 
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Figure 11. Heat Transfer and Flow Friction Characteristics of Solar 1. 
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Figure 12. Heat Transfer and Flow Friction Characteristics of Solar 4. 
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Tigure 13. Heat Transfer and Flow Friction Characteristics of Solar 6. 
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cure 14. Heat Transfer and Flow Friction Characteristics of 
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Figure 15. Heat Transfer and Flow Friction Characteristics of 


Brass Core, L= 0.5" 
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Figure 16. Heat Transfer and Flow Friction Characteristics of 
Brass Core, L=0.75". 
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Figure 17. Heat Transfer and Flow Friction Characteristics of 
Brass Core, L =1.0" 
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Figure 18. Heat Transfer and Flow Friction Characteristics of 
srass Core, L =1.5" 
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Figure 19. Heat Transfer and Flow Friction Characteristics of 
Brass Core, L = 2.25" 
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Figure 20. Heat Transfer and Flow Frict 


Brass Core, L= 3.0* 
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Figure 21. Combined Heat Transfer Characteristics of all Brass Cores. 
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Figure 24. Comparison of Brush and Moseley Recorders. 
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TABLE XIT 


SUMMARY OF GEOMETRICAL AND PHYSICAL PROPERTIES 


Matrix Material 


Specific Heat, Cy 


Thermal Conductivity, K. 


Material Thickness 


Hydraulic Diameter , D, 


Brass 0.5" 


Flow Length, L 


Frontal Area, A 


Volume, V 


Free Flow Area, A 
Conduction Area, A. 


Heat Transfer Area, A 


tr 


c 


Matrix Density, ? 


Porosity, p 


Compactness, G6 


Weight, W, 


Brass 0.75" 


Flow Length, L 


Frontal Area, A 


Volume,V 


9 


Free Flow Area, A. 


Conduction Area, A. 


Heat Transfer Area, A 


BRASS L/D, TEST CORES 


m 


Constants Common to all Cores: 


(Btu/lbm deg F) 


(Btu/hr ft deg F) 


(inches) 


(feet) 


(feet) 

(sq ft) 
(cu ft) 
(sq ft) 
(sq ft) 
(sq ft) 


(lbm/cu ft) 


(sq ft/cu ft) 


(1bm) 
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70-30 Brass 
0.092 
57.0 
0.001 


0.0013998 


0.04167 
0.05875 
0.002447 
0.051698 
0.007052 
6.9650 
61.7743 
0.87996 
2846.34 


0.15302 


0.0625 
0.05738 
0.003586 
0.050537 
0.006842 


10.2037 


Matrix Density, p m 
Porosity, p 
Compactness, (9 
Weight, W 
Brass 1.0" 

Flow Length, L 
Frontal Area, Ag, 
Volume, V 

Free Flow Area, A. 
Conduction Area, A. 
Heat Transfer Area, A 
Matrix Density, p m 
Porosity, p 
Compactuess, § 


Weight, me 


Brass 1.5" See Figure 10. 


Brass 2.25" 
Flow Length, L 
Frontal Area, Ag, 
Volume, V 

Free Flow Area, A. 
Conduction Area, A. 
Heat Transfer Area, A 
Matrix Density, Pm 
Porosity, p 


Compactness, § 


Weight, We 


TABLE XII (Cont) 


(lbm/cu ft) 


(sq ft/cu ft) 


(1bm) 


(feet) 

(sq ft) 
(cu ft) 
(sq ft) 
(sq ft) 
(sq ft) 


(lbm/cu ft) 


(sq ft/eu ft) 


(1.bm) 


(feet) 

(sq ft) 
(eu ft) 
(sq ft) 
(sq ft) 
(sq ft) 


(lbm/cu ft) 


(sq ft/cu ft) 


(1bm) 


7? 


S770 
0.8807 
2845.43 


0.20776 


0.08333 
0.060833 
0.005069 
0.05378 
0.007052 
i 5 
59a2oeo 
0.884076 
2748 .08 


0.30035 


0.18750 
0.060625 
0.011367 
0.05357 
0.007052 
31.3425 
64.8924 
0.883628 
LT eae 


Os7570 


Brasis 3.0" 
Flow Length, L 
Frontal Area, A, 
iG 
Volume, V 
Free Flow Area, A 
fr 
Conduction Area, A, 
Heat Transfer Area, A 
Matrix Density, Pm 
Porosity, p 
Compactness, B 


Weight, - 


TABLE XII (Contd) 


(feet) 

(sq ft) 
(cu ft) 
(sq ft) 
(sq ft) 
(sq ft) 


(lbm/cu ft) 


(sq ft/cu ft) 


(lbm) 
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O22 50 
0.060625 
0.0151516 
0.053573 
0.007052 
41.79 
60.5958 
0.883628 
Zhe 2 


0.9184 


L(inches) 


0.5 


1.0 
r.5 
Bao 


3.0 


TABLE XIII 


SUMMARY OF SLOPES OF BRASS TEST CORE £ AND j DATA 


WITH CORRESPONDING L/D, 


L/D 


Bos 
44. 
a. 
Ge . 
134. 


r/7@. 
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Slope (f) 


-0. 
-Q. 
-Q. 


-0. 


885 
o> 
or 


965 


941 


. 900 


Slope (j) 


-0.45 

-0.556 
-0.607 
07832 
-0.802 


-0.796 


APPENDIX A 

Description of Equipment 

The equipment used for the transient testing of core samples con- 
sists of a flow straightening device and air heating system which pre- 
cedes the matrix test section. Pressure taps measure both static pres- 
sure and pressure drop for the test section and the orifice for flow 
measurement. Thermocouples measure the temperature response of the 
matrix exit and the air temperature at the orifice. An ASME standard 
orifice is used to measure the flow rate, and a prime mover provides 
the air supply. 
Air Supply 

Air is used as the working fluid and is drawn through the equipment 
by a 30HP, multistage, Spencer Turbo-Compressor, which is rated at 550 
cfm using 220 V a.c. 
Flow Metering System 

An ASME standard orifice section using D, D/2 taps and thin edged 
concentric orifices was used for flow measurement. Orifice diameters of 
0.775 in., 1.232 in., and 1.540 in. were used with a 3.08 in. diameter 
pipe. 
Air Heater System 

The air heaters are bakelite frames wound with .0031 in. diameter 
nichrome wire spaced 1/32 in. apart, 50 to 52 wires per heater. There 
are two heaters per frame and fourteen frames totaling 28 heaters, which 
are wired in parallel across the input voltage source. Each pair of 
heaters is controlled by an individual selector switch so that depending 
on the flow rate, sufficient heaters may be used to give a twenty degree 


temperature rise to the incoming air. 
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Matrix Holder and Test Section 

The matrix holder and test section are constructed of closely 
machined polyethylene plastic; the tight fit minimizing air leakage. 
This section consists of a casing and a sliding drawer to hold the matrix 
under test. The casing has a removable frame on which the upstream thermo- 


couples, t,, are mounted and the upstream and downstream static pressure 


2 
taps. The sliding drawer contains the matrix and a plastic movable frame 


with the downstream thermocouples, t The flow channel is 3-1/16" by 


3° 
3-1/16" and can hold matrices of up to three inches in length. Matrices 
were placed in the holder and were surrounded on all four sides by styro- 
foam plastic insulation. This insured a snug fit of the matrices and 
also lessened heat loss from the matrix to the holder. 
Inlet Cone and Flow Straightener 
This section was designed by Piersall (13) and provided a uniform 
velocity profile to the air entering the heater section. 
Pressure Measuring System 
Pressure taps are located in the matrix holder upstream and down- 
stream of the matrix and in the pipe at one diameter and one half dia- 
meters on either side of the orifice. Each pressure tap is connected 
by flexible tubing to its corresponding manometer and draft gage. The 
following instruments were used: 
1. Ellison Draft Gage Company, 0-3" inclined gage. 
2. Ellison Draft Gage Company, 10" manometer. 
3. Ellison Draft Gage Company, 20'' manometer. 
4. Merriman Instrument Company, 120" manometer. 


5. Precision Thermometer and Instrument Co., mercury barometer. 
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Temperature Measuring System 


Temperatures are measured in the system at the orifice, at the 
inlet to the system, between the heaters and the matrix and at the matrix 
outlet. Thermocouples are used for all measurements. The thermocouple 
at the orifice is a single 30 gage thermocouple, referenced to an ice 
water junction (i.e., 32 deg F), and output read on a Rubicon portable 
potentiometer. 

There are four sets of 30 gage iron-constantan thermocouples, each 
set consisting of five thermocouples connected in series. These thermo- 
piles were constructed and wired by Traister (17). Two of the four sets 
of thermocouples had each thermocouple individually wrapped with teflon 
tape to prevent shorting and were placed in an aluminum tube mounted in 
a frame at the exit of the inlet cone. The aluminum tube acted as a 
radiation shield to prevent the heaters immediately downstream from affect- 
ing the readings. Each of these sets were designated th and measured 
the temperature of the incoming air. The third set, designated to» was 
placed in a frame between the heaters and the matrix. This set was 


"bucked" against one of the t, sets so that the output of the two sets 


1 


measured the difference between t, and t.- In a similar manner, the last 


thermopile was wired to the second set of ty and measured the dif- 


ference between incoming air and the matrix outlet, ty = tl: For data 


taking, the output of t, - t) was led to one channel of a Hewlett- 
Packard, Moseley Division, Model 7100B dual channel strip chart recorder. 
This thermocouple output was also led into an Astrodata model 886 Wideband 
Differential D.C. Amplifier where it was amplified 1000:1. 

Differentiator 


The amplified thermocouple output from the D.C. amplifier was fed 


to the differentiating circuit on the analog computer. The computer used 
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was a Donner Model 3500 portable analog computer. The circuit used is 
shown in the theory section. The differentiated thermocouple response 
was then fed into the remaining channel of the strip chart recorder where 
the derivative could be compared directly with the undifferentiated signal. 
Heater Power 

Power for the heaters was supplied from the 250 V D.C. source in 
the laboratory. For low flow rates, the voltage at the supply panel was 


adjusted to give the proper temperature rise to the air. 


85 


APPENDIX B 

Data Reduction Relationships 

This appendix summarizes those data reduction relationships of 
importance in calculating flow friction and heat transfer values from 
the collected data. 
Geometry 

Accurate determination of dimensions and physical constants are 
necessary to minimize error so that comparisons between cores may be 
made. Compact heat transfer surfaces use three geometric parameters 
that allow comparison to be made between matrixes. These are: 

1. Hydraulic Diameter 


2 _ 4 x free flow area _ 
a oe heat transfer area oo a 


2. Porosity 


free flow area 


frontal area 7 A./Ae, (B-2) 
3. Area Compactness 
. heat transfer surface area _ 
p= matrix volume A/(A, 1) Cae 


Dividing (B-2) by (B-3) results in: 


r= PB (B-4) 
Through this equation, knowing any two of the above relations leads 


to the third. 


Maximum Slope 


The following sketch of the generalized cooling curve will serve to 


help explain the maximum slope relationships: 
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Given that the maximum slope of the above curve is a function of 


N., and A. zy A ait 
| = bm) 
i ) 
where: ea 


T = generalized time variable ~ hA@/(W_C.) 


N = hA/(mc_) 
tu EP 


Yn, me ,/(W.c,) x @ 


and: 
of (at )= ™m Cp dg (B-5) 
Neu Ws Cs 
Furthermore 


Cr. i tT; en th - Lf, 
fy, - ti = rc - tz Ae 


and its derivative is: 





%-t; / 
a = ol 27 os (B-6) 
Ce, = ‘) (ty 8) (% “) 


Combining equations (B-5) and (B-6) gives us: 
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eet -) a Was | d (t,.-t 








MAX 


From the generalized curve: 


\2 (tee 2) e 


x/(chart speed) = dO 
(too so" e1? . 
Sey '? 
combining with: 
: _ matrix capacity a -l 
eaten flow stream capacity rate : IC a 


and equation (B-7): 


d ae Ca C¢ ) 
Ate = til = fas = (Y)chare speed 
a io ) Bs (B-8) 
Mtu MAX 


This value of maximum slope and A are then used to enter Table I 
or Figure 1 to get the corresponding value of Ney 
Maximum slope from analog output: 
The analog computer reads dy/dx directly in millivolts per second. 
This value, along with (m) in millivolts produces maximum slope from 


the following equation: 
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ul ee =| 





te - & 
= meee ‘ = 5 . dy/dx (B-9) 
ol me ) Me), m 
Neu MAX 


Flow rate 
The mass flow rate was calculated from ASME Power Test Code (1) 


as modified by Murdock (13) by the following equation: 


a 350K a* F el Ps (B-10) 
0 a Oo 


where 


Ss 
Vi- BF 


aA 
il 





C = orifice coefficient of discharge (12) 


§S = ratio of orifice diameter to pipe diameter 


d = orifice diameter in inches 


Fs = orifice plate thermal expansion factor 
= fluid thermal expansion factor 
= specific weight of fluid flowing = P/(RT) assuming 
perfect gas. 
P = absolute static pressure at orifice (lbf/sq ft) 


R = Gas constant for air: 53.35 (ft-1b,/1b_ degR) 


T = absolute temperature at orifice(deg Rankine) 


re = pressure drop across the orifice in inches H,0 


Substituting the expressions for K and qT in equation (B-10) yields: 





Mm =359=8 do Fa Aip, © (B-11) 


1-4 RT 


89 


From (1), Fig 40A 


\- (0.4) + 0.358") 


YX 


k = 174? for aifjeatio of fe sc 
Ce v 


_ AP, (in Hg abs) 
P (in Hg abs) 


also from (1) fig. 38 


F = 1.0 
; P 
S 
P= (Patm - jac) (04912 */44) IbF/ p42 
P sen = local atmospheric pressure in inches Hg 
P = static pressure upstream of the orifice plate in inches Hg 


O 


Substituting into (B-11) the above expressions with the necessary physical 


constants to make the equation dimensionally consistent yields: 








m= 589.0) —_Ce @fi-(41+0354)| 4% | _1 
V/-3t | fn) 99.02 
13.6 
\ B-12 
AR ( Pat m ~ 8, “ ; 
t, + 459.7 


Reynolds Number 


Reynolds Number is defined as: 


A nee (B-13) 


where G is the mass flow velocity based on the free flow area, A. 
G = m/A_ = m/(pA B-14 
m/A. = m/(p er ( ) 


substituting: 


and from (B-4), r= P/B agai Di = ay 
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therefore: 


N= a (B-15) 
R MAB 
Substituting (B-3) into (B-15) 4s A 
yielded. epvel 
= (4mL)/(A 4) (B-16) 


Fanning Friction Factor 
The following equation and the accompanying sketch describes the 
flow system under consideration: 


Rs ‘ A 2: [(Kc rine 2(G-1) + FB Yon (1-6 key 


ENTRANCE Flow ea Exir “amen 
EFFECT ACCELERATION FRICTION fa FFECT 


LLIDM MILL AMAA LLALMLH ALLL ALMA MMMMMML AMAL ALAA MAAYG SL | | 


i 
' ' 
b 
; 
j LALLILLULLILMLLIA LIM MMALLA LEAMA ELIMI TTA L TEA | 
| VLA TLLLETET ETAL ALLTEL TIL TTT ALT UILLTA I 
| i 
FLOW CLE T TTT AT TLTTL TLL THLA TILA ATLL TUTTI ETD 
{ { 
| 
tt TTT TT TEE 


ia oS 2 


Because a gas is the working fluid the changes in pressure from 1 
to a and from b to 2, in the preceding sketch, are small in comparison 


with the total pressure; consequently, u,= 4, and us Uy: Since 


oe 


testing is done at moderate temperatures and pressures, the perfect gas 
law is assumed valid. From (B-l), A/A = L/r, andv=1/e ; substitut- 
ing into equation (B-17) and solving for f yields, for the isothermal 


f=|2g- Pm (22) -(Se + Ke) Refs — Hehe (7 4 -£)ir 9) 2 (B-18) 


where Gy = G20 and subscript 1 refers to upstream values and 
subscript 2 refers to the downstream values. KW and K, are the entrance 
and exit loss coefficients respectively and are dependent on porosity, 
geometry of flow cross-section and Reynolds No. within the core (8). 
These coefficients were calculated from the analytical expressions derived 
by Kays (8). 
C =Contraction Ratio = 0.611 + .045p + .344p mie 
Ka = Velocity distribution factor 
1. Circular tubes 
laminar flow 
= 1.365 

turbulent flow 

f = 0.049N (Fanning friction factor, circular tubes) 

Ky = 1.09068(4f) + 0.0588 4£ +1. (B-19) 


2. Gap, laminar and turbulent flow 


ae =0.6 (B-20) 


Ky (circular) - l 


3. Square, laminar and turbulent flow 


K , (square) - l 


K (circular) — hed? (B-21) 
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4, Triangle, laminar and turbulent flow 


K (triangle) - l 
aS CO 128 (B-22) 
Ky (circular) - l 


Using the appropriate K, for the passage cross-section and the 


d 
appropriate flow, 


2 
Oe ee (B-23) 
l - 2c + cs (2K, - Ll) (B-24) 
Ko 2 
C 
Cc 


Using an order of magnitude approximation, for small pressure differ- 
entials the first term in (B-18) is the greatest contributor to the 


friction factor. The approximation: 


a : =) ae Pa P= P, reduces (B-18) to: 
f= [2ge A, 42 ~ (ker Ke) - 22 (inet) TH a2) 
G* Pam ‘a 


By substituting (B-1l) and (B-14) into (B-25) it can be shown that: 
£ >/g 
Colburn j-Factor 
The Colburn j-factor is defined as: 


- Ee ee. 2/2 
J = Nge “pr Ge ‘Pr ay 


Substituting (B-14) for G and multiplying by A/A yields: 





A 
= hA S 2/3 > 2 oe 
j fie - 7: N,,, - but Neu hA/(me |) 


therefore: 


> Ze 
j =N a (A /A)N, . 
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and by substituting equations (B-1) and (B-4) it can be shown that 
j P/B 

Heat Transfer Power and Friction Power 

Relative performance of matrices under comparison can be determined 
by an evaluation of the heat transfer power vs. flow friction power. 

The higher the plot of hee vs. ee the better the core (9). 

The heat transfer power per unit area per degree temperature differ- 
ence is: 

h = ST als (4) Np J (B-27) 


Evaluating er , and Noy at standard conditions of 


500 deg F and one atmosphere for convection (9), 


0.2477 Btu/1lb deg F 


0.0678 lb/hr ft 


o 
p 


Noy 


Equation (B-27) evaluated at standard conditions becomes: 


0.0413 lb/ft? 


0.671 


- 0.07195 (44 (Ne J) re Of (B-28) 


The flow friction power per unit area is (9): 


«= flag 23. tL. rT Ne | a 


Evaluating equation (B-29) at the Standard conditions shown above gives 


the equation for E ota’ 


-7 A rr ey 
E etd = 1.11 x10 (é. 1000 ae, (B-30) 


94 


APPENDIX C 
Digital Computer Program for Data Reduction 

The digital computer program used by Traister (16) was written for 
a CDC model 1604 computer. The acquisition of an IBM 360 Systen and 
attendant removal of the CDC 1604 computer required the program to be 
converted from Fortran 60 language to Fortran IV language. This program 
takes the sample core geometry and all raw data and calculates all the 
heat transfer and friction résults used in this report. 

The program uses a curve fitting interpolation subroutine to deter- 
mine the value of New from maximum slope and conduction parameter. This 
subroutine uses Howard's conduction parameter data. (Table I). While 
not used in this investigation, a cyclic method for finding the value 


of Neu exists and a subroutine for this method is included in the program. 
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